Ortho— H2 and the age of interstellar dark clouds 



Laurent Pagani 

LERMA, UMR8112 du CNRS, Observatoire de Paris, 61, Av. de Wbservatoire, 75014 Paris, France. 

laurent .pagani@obspm.fr 
and 

Evelyne Roueff 

LUTh, UMR8102 du CNRS, Observatoire de Paris, 5, place Jules Janssen, 92195 Meudon cedex, France 

and 

Pierre Lesaffre 

LERMA, UMR8112 du CNRS, Ecole Normale Superieure, 24 rue Lhomond, 75231 Paris Cedex 05, France 



ABSTRACT 

Interstellar dark clouds are the sites of star formation. Their main component, dihydrogen 
exists under two states, ortho and para. H2 is supposed to form in the ortho:para ratio (OPR) 
of 3:1 and to subsequently decay to almost pure para-H2 (OPR < 0.001). Only if the H2 OPR 
is low enough, will deuteration enrichment, as observed in the cores of these clouds, be efhcient. 
The second condition for strong deuteration enrichment is the local disappearance of CO, which 
freezes out onto grains in the core formation process. We show that this latter condition does not 
apply to DCO+, which, therefore, should be present all over the cloud. We find that an OPR > 
0.1 is necessary to prevent DCO+ large-scale apparition. We conclude that the inevitable decay 
of ortho-H2 sets an upper limit of ^^6 million years to the age of starless molecular clouds under 
usual conditions. 

Subject headings: ISM: abundances — ISM: clouds — evolution — ISM: molecules — astrochemistry 



1. Introduction 

The low mass star formation process is rela- 
tively well understood, however its first step de- 
tails remain uncertain, especially regarding the 
question of time and age. Still much debated 
questions are: how long it takes for a cloud to 
form from the diffuse atomic medium, how long 
a molecular cloud lives, how long it takes to form 
a condensation that will evolve into a prestel- 
lar core that subsequently collapses to form a 
protostar. Today, much debate exists upon the 
lifetime of clouds, e.g. iHartmann et al.l ( 2001 ) 
defending a short lifetime (1 to a few m illion 
years) while Tassis fc MouschoviasI ( 20041) and 
Mouschovias et al. ( 20061) claim a tvpical age of 



10 million years, but most arguments are either of 
limited statistical significance or model-dependent 
(unknown magnetic field strength, small scale tur- 
bulence rapid dissipation, etc.). Clues are needed 
and while the usual chemical age modeling is not 
satisfying (problem of unknown initial conditions) , 
we present here a simple constraint based on basic 
chemistry: the absence of DCO+ in dark cloud 
envelopes can only be explained if the clouds are 
young enough and we propose to determine an up- 
per limit for this age which is as much independent 
from the initial conditions as possible. 
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2. The abundance of CO and DCO+, and 
the deuterium chemistry 

Dark clouds contain dust grains embedded in 
interstellar gas, itself mostly molecular as these 
clouds are self-shielded against our Galaxy's UV 
background. These dark clouds ar e cold (typi- 
cally 10 K, with cores as low as 7 K, Pagani et al.l 
2003 . 20051) . and deep inside, where the density 
reaches a few 10* cm~'^, the heavy species freeze 
out onto grains to form ices. Both these ices 
and the gas are subject to a rich chemistry. Re- 
cently, strong deuteration (deuterium enhance- 
ment with respect to hydrogen carriers: DCO^ 
versus HCO+, HDCO and D2CO versus H2CO, 
etc.) has been recognized as a useful tool to 
better understand both t he chemistry itself and 
the star formation pro cess ([Bergin fc Tafallall2007t 



Ceccarein et~al]l2007t ) 



H2 is the main hydrogen reservoir in dark 
clouds, and similarly, HD is the main deuterium 
reserv oir. Its relative abundance to Ho is ^3 x 
10-5 ()Hebrardll2006t iLinskv et al.lbood) . IWatson 



(|1976[ ) was the first to invoke chemical deuterium 
enrichment via the important reaction 



H+ + HD ^ H2D+ + H2 + AE 



(1) 



AE — 232 K when all species (reactants and 
products) are in their ground state. In cold dark 
clouds, the forward reaction is strongly favorable, 
creating H2D+/H^ ratios much larger than the 
original HD/H2 ratio. If CO is abundant (X[CO] 
« 1.5 X lO"'*, where X[] denotes the relative abun- 
dance of a species with respect to H2), H;^ reacts 
5 times more often with CO than with HD, and 
forms HCO'*', strongly decreasing the abundance 

of H2D+. 

Once H2D+ is formed, it can be enriched in deu- 
terium by reacting with HD again to form D^H^ 
and Vlit, as first recognized by lRoberte et all (|2003h 



HoD^ 



HD 



D2H+ + H2 + 187K 



(2) 



D2H+ + HD ^ D+ + H2 + 234K (3) 

These two reactions are again favorable in the 
forward direction at low temperatures because of 
their exothermicity. Then, these Hj]' isotopologues 



can easily transfer their deuterons to other species: 

CO + H2D+ ^ ^HC0+ + iDC0+ + iH2 + ^HD (4) 

CO + D2H+ ^ iHCO+ + ^DCO+ + iHD + ^D2 (5) 

CO + D+ ^ DCO+ + D2 (6) 
1 



N, 



H2D+ 



N2 + D2H+ 



-N,H^ 



:N,D+ + iH2 + ^HD (7) 



rNoH^ 



3 

f ^N2D+ + IhD + ^D2 (8) 
N2 + D+ ^ N2D+ + D2 (9) 



When CO is abundant, H2D+ is much less abun- 
dant as noted above and also reacts with CO 
preferentially, quenching all the other deuteration 
paths. CO also reacts with species like N2H+ and 
N2D"'" to destroy them: 



N2H+ + CO ^ HCO^ 
N2D+ + C0 ^ DCO^ 



N2 
N2 



(10) 
(11) 



The strong correlation betwee n CO freeze-out onto 
grains and deuteration (e.g. ICrapsi et al. 2005) 



confi rm the chemical models (e.g. [Roberts et al 



20031) as sketched above and has led to the conclu- 
sion that CO freeze-out was necessary for strong 
deuteration to occur. There is however one ex- 
ception to this statement: the deuteration of 
CO itself can form DCO+. In first approxima- 
tion, if H;^ is mainly destroyed by CO, X[H2D+] 
varies like XfCO]""'^, whereas the production of 
DCO+ is proportional to X[CO] and X[H2D+] 
(reaction 4). Thus the dependency upon X[CO] 
cancels out and X[DCO+] should remain approx- 
imately constant. Figure [1] shows the abundance 
of X[DCO+], X[H+ isotopologues] and X[N2D+] 
as a function of X[CO] in a steady-state chemi- 
cal model. The X[Hi^ isotopologues] is the sum 
of the abundances of all 3 deuterated Hj]" iso- 
topologues, which represent altogether the main 
deuteration partners of spec ies like CO and N2 . 
The model is the same as in IPag ani et al. I iiooi) 
but CO now runs from an undepleted relative 
abundance of 1.5 x lO^'' down to a depletion by 
a factor 300 (X[CO] = 5 x 10"^), for a cloud 
density of lO'* cm'^ at 10 K. When X[CO] starts 
to drop, X[H3 isotopologues] increases faster than 
proportional (indicated by the -1 slope) because 
the D2H+ abundance depends both on the de- 
struction of Ht by CO and on the destruction of 
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H2D+ by CO, and the abundance depends 
also on the destruction of D2H+ by CO. This ex- 
plains that at the beginning, X[DCO+] increases: 
the deuteration capabilities increase faster than 
X[CO] decreases. While X[CO] keeps decreas- 
ing, the deuteration capabilities reach saturation 
and X[DCO+] becomes proportional to X[CO] 
(traced by the -1-1 slope). X[N2D+], a represen- 
tative gas-phase deuterated species, follows the 
deuteration capabilities of the cloud, as expected. 
It is clear that DCO+ has a behavior different 
from the other deuterated species, with an abun- 
dance varying only by 50% while the CO abun- 
dance drops by a factor of 30, and most impor- 
tantly by reaching its maximum abundance when 
CO is high and deuteration capabilities are down. 
Therefore, one would expect DCO"*" to be present 
everywhere in the cloud. It is not the case. Indeed, 
DCO+ is, as the other deuterated species, only de- 
tected in the de pleted prestellar cores on extents 
simil ar to NH3 ( Butner et al. 19951 lJuvela et al. 
20021 ) and nowhere else. This can only be ex- 
plained if a high II2 OPR is present. Indeed the 
II2 OPR is an important controller of the cold 
clo ud deuteration chemistry , as alr e ady discussed 
in iPineau des Forets et al.l (1991): iFlower et al 



(|2006l ) and lPagani et all (Il992ll2009[) . 



3. The role of ortho H2 

II2 has 2 different possible nuclear spin states 
due to the proton spin of ^ : parallel (ortho, I 
= 1, weight 21+1 = 3 ) and anti-parallel (para, 
1 = 0, weight = 1). The ortho state corresponds 
to the odd rotational levels (J — 1, 3, 5, ), the 
para state to the even levels. The first ortho state 
(J = 1) is 170 K above the ground (J = 0) para 
state. The ortho-para exchange is possible neither 
via radiative processes, nor via inelastic collisions. 
It can occur in the gas phase and on the surface 
of dust grains (Le Bourlot 2000,). The conversion 
rate on a solid surface i s however highly unce rt ain 
(jWatanabe et al.lbOloHchehrouri et alj|201ll) and 



it is most probable that this is mainly realized 
by g as-phase chemical rea, c tions with H"*" or Hj" 
ions toalgarno et al.' ll973l: IPlower fc Wattlll984l 
ILeBourlot 1991; FlowerelalJl2006|). Theimpor- 
tance of ortho-H 2 in the deuteration pro c ess wa s 
first studied bv IPineau des Forets et al. ( 1991 ): 
Pagani et al. ( 19921 ). and later on. Flower et al 



regulator. However, they considered the conse- 
quences of the OPR evolution only for the prestel- 
lar core format ion itself. In a, stud y of LI 83 main 
prestellar core, jPagani et al. I (I2OO9I) noted that the 
H2 OPR evolution timescale required by the ob- 
served N2D+/N2H"'" ratio was on the order of the 
free-fall time. They also noted that the OPR had 
to drop below 0.01 for N2D"'" to become detectable. 

The II2 OPR value is not easily accessible. H2 
does not emit in cold dark clouds and is seen only 
in shocked regions, especially protostar outflows, 
or in absorption in UV in low extinction regions. 
The energy release from the H2 formation reac- 
tion, assumed to be 1.5 eV in the equipartition 
hypothesis, is large enough to populate many lev- 
els of both ortho and para species: it is hence 
generally accepted, though not certain, that the 
OPR is 3:1 when II2 forms on grains. Once re- 
leased in the gas phase, II2 tends to relax towards 
its thermal equilibrium state but never reaches it 
because at trace abundance levels, ortho-Il2 de- 
struction is compensated by fresh H2 formed onto 
grains from residual H and from destruction of 
ions such as Ht and Il2D''~ that can release an 



ortho-H2 (|Le Bourlotlll99ll : lFlower et al.ll2006l ). It 

is therefore necessary to model the chemical evo- 
lution of the II2 OPR to understand the evolution 
of dark clouds. Once II2 is formed, it can react 



with H"*" and to exc hange states (jLe Bourlot 



1991: Flower et al,l 



. For example. 



0-H2 + 



P H2 



AE 



(12) 



(|2006l ) revealed its importance as a deuteration 



AE — 170 K is the energy released as II2 goes 
from J=l to J=0 levels. The reverse endothermic 
reaction is thus difficult to obtain in a cloud at 
10 K and ortho-H2 slowly decreases to low abun- 
dance levels (wl0~^ - 10"''). Similarly, in reac- 
tion 1, except for HD, all the species have ortho 
and para states and AE therefore depends in fact 
on the considered initial and final states. If only 
para states are involved in the products, AE is 
maximal (AE = 232 K with para- H+ and 265 K 
with ortho-H;^, HD being always in its ground 
state, J = 0) but if both products are in their or- 
tho ground state while the reactants are in their 
lowest state, para-Hg" ground state and HD (J 
— 0), AE becomes negative and the reaction is 
slightly endothermic (H2D+ lowest ortho state is 
87 K above the para ground state, and the total 
required energy is 25 K). Similarly, while the re- 
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verse reaction is impossible in a 10 K cloud with 
both para reactants, it becomes rapid with ortho- 
H2D+ and ortho-H2 even at lOK, as this chan- 
nel is slightly exothermic (AE — 25 K) to pro- 
duce para-H;^. A few such channels are opened 
with ortho-H2 allowing the efficient destruction 
of the deuterated isotopologues. Ortho-H2 
is thus an important chemical poison to deuter- 
ation in dark clouds. Its abundance is therefore 
critical to control the deuteration en hancement 
([Flower et al.ll2006t [Pagani et al.ll2009( ). Because 



000+ is not observed everywhere in the clouds 
outside the depleted cores where it i s detected 
(jButner et all Il995l : ljuvela et all l2002[ ) while its 



abundance should be maximal, this sets a mini- 
mum level for the abundance of ortho-H2. In turn, 
this minimum level gives rise to an upper limit to 
the age of the cloud, as ortho-H2 abundance must 
eventually decrease to low levels as implied by re- 
action (12), and in order to explain the deuteration 
enhancements seen in the depleted cores, which. 



tion up to 10^^, Li 


s et al. 2002 


2002: Parise et all 


2004h. 



2002; van der Tak et al 



The chemical model 



To address these questions, we used the deuter- 
ation network from iRoueff et all (^005., .2007.) in 
which we have included the ortho-para chemistry 
' I (l2009l) . based itself on 



as described in Pagani et al. 



the w ork by ,Walmslev et al.l (|2004l ) ; iFlower et al 



( 2006 ) and the new rate coefficients of lHugo et al 



20091 ). We have also included analytical approx- 



imations of the spin state-dependent dissociative 
recombin ation rates of H j" iso topologues from the 



tables in Pagani et al. 



(|2009[) . The deuteration 
network includes all simple species (up to 6 atoms) 
based on C, N, O, and S with up to 5 deuterium 
substitutions (CD^). For those reactions with a 
small endothermicity (E/k^ < 600 K), we have 
corrected the endothermicity when the reaction 
with ortho-Il2 was involved (e.g., CD -I- II2 — > CH 
-|- HD is endothermic with para-H2 and exother- 
mic with ortho-112). For all reactions with rare 
species producing II2 (e.g., HCN+ + CN+ + 
II2), we have always assumed that II2 is released 
in its para form, the most favorable energetically 
(similarly, only ortho~D2 is considered in similar 
reactions with deuterated species). This does not 
noticeably accelerate the decay of ortho-H2. We 



checked that after, e.g., 10^ years, the para-Il2 
production rate is 99.9% dominated by reactions 
of ortho-Il2 with H;^ and 11+ . Our model is not 
aimed at reproducing accurately a given cloud but 
it allows to follow the evolution of the II2 OPR 
in a molecule rich gas. We have computed the 
chemical evolution of a cloud of constant density 
(n(Il2) = IC'cm"^), constant temperature (10 K), 
average grain radius of 0.1 /im and cosmic ion- 
ization rate of 1 x 10~^^s~^, which represent the 
usual conditions met in starless dark clouds. We 
have started with all species being atomic except 
H which is considered to be already entirely con- 
verted to II2 and D to HD and we have consid- 
ered different OPR starting ratios from 3 down to 
3 X 10"'^ to account for the possibility of H2 pro- 
duction OPR below 3. Figure 2 shows the evo- 
lution of DCO"*" and ortho-Il2 with time. We 
have traced the detectability limit of DCO+ for 
a cloud with a H2 column density of lO^^cm"^, 
considering that a DCO+ J:l— >0 line with a width 
of 0.5 km s^^ and an intensity of O.IK would 
be easily detectable with present day radiotele- 
scopes. With the RADEX non-LTE radiative 
transfer model ( van der Tak et al"1l2007l ). we find 
that a DCO"*" column density of ~2 x 10^^ cm~^ 
is detectable, equivalent to a relative abundance 
to H2 of ~2 x 10~^^. Figure 2 shows that for 
any initial H2 OPR below 0.1, DCO+ should be 
detected throughout the cloud in less than 3 x 
10^ years. For an initial OPR = 0.1, the DCO+ 
becomes marginally detectable around ~5 x 10^ 
years, disappears and finally becomes clearly de- 
tectable after 2 My. For OPR > 0.1, it takes 3 to 
6 million years to become detectable, i.e. when 
the H2 OPR drops below '--^0.03. Compared to 
the H2 OPR drop from 3 to ~0.01 in depleted 
cores which takes le ss than 2 x 10^ years to happen 
( Pagani et al. 20091 ). the difference is twofold: as 
density increases during the prestellar core forma- 
tion, 1) the chemistry accelerates, and 2) depletion 
sets in and H+ is no more destroyed by the other 
species which have now disappeared from the gas 
phase (H;^ becomes also more abundant but to a 
lesser extent). H+ can even become the dominant 
ion. The higher abundance of H+ and H'^ and the 
accelerated chemistry provoke the rapid decline of 
ortho-H2 . 

On the same figure, we also trace the abundance 
of the sum of the deuterated isotopologues in 
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the case H2 OPR = 3. It shows that under 1 My, 
deuteration is dominated by the carbon chemistry 
(CH2D+, C2HD+, etc.) which is not dependent 
upon the ortho-H2 abundance due to exothermic- 
ities close to 400 K. It is only beyond 1 My that the 
H;^ deuteration network takes over when ortho-Il2 
starts to be negligible. If ortho-H2 is low since the 
beginning, then the Hj^ isotopologue contribution 
to the abundance of DCO^ is important from the 
start which explains the disappearance of the dip 
between 1 and ~3 My. 

The model is also sensitive to two other ini- 
tial conditions: the abundance of metals and the 
cosmic ray ionization rate. Metals, especially al- 
kali and alkaline earth metals like Ca, Na, K, and 
others like Fe, have a direct influence on the elec- 
tronic equilibrium of the chemical model which in 
turn changes the dissociative recombination effi- 
ciency for DCO'*' and therefore its abundance. For 
standard C (1 x 10"" s'^) and H2 OPR (3), we 
have varied their abundance from 3.4 x 10~® to 
1.3 X 10^^. For the lowest abundance (50% of 
the default one), DCO+ becomes detectable much 
earlier ( 2 x lO'^ y) than in the standard case, 
drops below the detectability limit after 1 My and 
is detectable again after 4 My. For the highest 
abundance (twice the default one), 000+ is de- 
tectable after 7.5 My. The cosmic ray ionization 
rate (Fig. |3]): if the rate increases, the produc- 
tion of 11+ and increases, therefore accelerat- 
ing the decay of ortho-H2 and increasing the abun- 
dance of HCO+ and DCO+. For C = 3 x lO'^^ 
s~^, DCO+ becomes detectable in less than 10^ 
years. Values of ( in between 1 and 3 x 10^^^ s^^ 
would bring our model in agree ment wit h some re- 
cent estimates of cloud l ifetime ([Enoc h et al. 200^ 
Hatchell fc Fulleill2008l) . Conversely, for C = 3 x 



10 



-18 



the chemistry is considerably slowed 



down, DCO+ becomes detectable only after 28 My 
and steady-state is reached only after 47 My. It 
must be noted however that such a low ioniza- 



tion r ate seems unlikely following iPadovani et al 
()2009l) . 



5. Conclusions 

Ortho-H2 abundance is most probably the 
main limiting parameter for DCO+ production. 
The absence of DCO"*" outside depleted cores, 
while its abundance is expected to be maximum 



for CO abundances close to 10~^, is due to a high 
ortho-Il2 abundance in dark clouds (OPR > 0.1) 
maintained long enough. Such high abundance 
can survive thanks to low and H^^ abundances, 
remaining low via charge or proton exchanges with 
many neutral species. However, ortho-H2 eventu- 
ally disappears. This sets an upper limit to the age 
of DCO^-quiet dark clouds of 3-6 million years 
after H2 has formed for the most probable ini- 
tial conditions. Higher cosmic ray ionization rates 
would increase the abundance of H+ and Hj and 
can thus only shorten the maximal possible age 
of the clouds to values which are too small while 
lower rates would tend to relax the age constraint 
but seem ve ry unlikely. Thi s resu lt strengthens 



the re s ults oflHartmann et al.l (120011) ; Enoch et al 



(|2008l) : iHatchell fc Fulleil (120081) and is no t com- 
pati ble with those of Tassis fc Mouschovias ( 20041 ) 



and Mouschovias et al 



(|2006[ ). From the chemical 
point of view, cold dark clouds seem to be really 
short-lived. 

We are grateful to Patrick Hennebelle for fruit- 
ful discussions. ER acknowledges the support of 
the ANR NO. 09-BLAN-020901. 
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Fig. 1.— DCO+, N2D+ and H;^ isotopologues 
abundances as a function of CO abundance in a 
steady state model. Slopes of +1 and -1 are traced 
to better visualize the variation slope of the differ- 
ent ions. DCO+, on the contrary to other deuter- 
ated species like N2D+ reaches its maximum abun- 
dance for high CO abundances. 
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Fig. 2.— Variation of DCO+ and ortho— H2 
abundances in a pseudo time-dependent chemical 
model for different starting H2 OPRs from 3 x 
10"'^ up to 3. The abundance of the sum of the 
deuterated H3 isotopologues is also traced as a 
thick gray line (for the case H2 OPR = 3). The 
horizontal dotted line in the lower box indicates 
the limit of detection of DCO"*" and the circles in 
the upper box mark the corresponding OPR val- 
ues. 
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Fig. 3.— Variation of DCO+ and ortho— H2 
abundances in a pseudo time-dependent chemi- 
cal model for different cosmic ray ionization rates 
from 1 X 10~^* up to 1 X 10~^^. The horizon- 
tal dotted line indicates the limit of detection of 
DCO'*', the arrows mark the time at which DCO"*" 
abundance crosses the detection limit. 
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